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Abstract 

In this paper I report results for simulations of the three-dimensional gauge glass and the four-dimensional XY 
spin glass using the parallel tempering Monte Carlo method at low temperatures for moderate sizes. The results 
are qualitatively consistent with earlier work on the three- and four- dimensional Edwards-Anderson Ising spin 
glass. I find evidence that large-scale excitations may cost only a finite amount of energy in the thermodynamic 
limit. The surface of these excitations is fractal, but I cannot rule out for the XY spin glass a scenario compatible 
with replica symmetry breaking where the surface of low-energy large-scale excitations is space filling. 
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There has been an ongoing controversy regard- 
ing the spin glass phase. There are two main the- 
ories: the "droplet picture" (DP) by Fisher and 
Huse [1] and the replica symmetry breaking picture 
(RSB) by Parisi [2]. While RSB follows the exact 
solution of the Sherrington-Kirkpatrick model and 
predicts that excitations involving a finite fraction 
of the spins cost a finite energy in the thermody- 
namic limit, the droplet picture states that a clus- 
ter of spins of size / costs an energy proportional 
to where > 0. It follows that in the thermody- 
namic limit, excitations that flip a finite cluster of 
spins cost an infinite energy. In addition, the DP 
states these excitations are fractal with a fractal 
dimension ds < d, where d is the space dimen- 
sion, whereas in RSB these excitations are space 
filling, i.e. ds — d. Krzakala and Martin, as well 
as Palassini and Young [3] (referred to as KMPY) 



find, on the basis of numerical results on small sys- 
tems with Ising symmetry, that an intermediate 
picture may be present: while the surface of large- 
scale excitations appears to be fractal, only a finite 
amount of energy is needed to excite them in the 
thermodynamic limit. 

The differences between DP and RSB can be 
quantified by studying the distributions P{q) and 
P(qi) of the spin overlap q and link overlap qi. For 
finite systems, the DP predicts two peaks at ±q-EA, 
where ^ea is the Edwards- Anderson order param- 
eter, as well as a tail down to g = that vanishes 
in the thermodynamic limit, like ^ for pertur- 
bations introduced by a change in boundary con- 
ditions!^. In contrast, RSB predicts a non-trivial 
distribution with a finite weight in the tail down 

^ In this work I use 8' instead of S as I introduce the 
excitations thermally with fixed boundary conditions [6] . 
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to g = 0, independent of system size. In addition, 
DP predicts the variance of the hnk overlap to fit 
a power law of the form 

Yniiqi) ^ a + bL^^'' (1) 

where a = and, as shown in Ref. [6], = 6' + 
2{d — ds), whereas for RSB one expects a > 0. 

While there has been considerable work on Ising- 
type spin glass systems, only few attempts [4,5] to 
understand the nature of the spin glass phase for 
models with a vector order parameter have been 
made. In this work I review recent results for the 
three-dimensional (3D) gauge glass and the four- 
dimensional (4D) XY spin glass. 

The Hamiltonian of the vector models analyzed 
can be summarized by 

(hi) 

where the sum ranges over nearest neighbors on a 
hypercubic lattice in d dimensions of size TV = L'^. 
Here 4>i represent the angles of the XY spins as 
S = (cos((/)), sin(0)) and therefore |S| — 1. Peri- 
odic boundary conditions are applied. In the case 
of the gauge glass, I set d = 3 and Jy = 1 for all 
The Aij are quenched random variables uni- 
formly distributed between [0, 27r] representing the 
line integral of the vector potential directed from 
site i to site j. For the 4D XY spin glass I set d = 4 
and Aij — for all The Jij are chosen accord- 
ing to a Gaussian distribution with zero mean and 
standard deviation unity. 

The spin overlap q is traditionally defined as 

1 ^ 

5-]^Esr-sf, (3) 

where a and /3 represent two replicas of the system 
with the same disorder. For this quantity to be a 
sensible order parameter, it has to be maximized 
with respect to all symmetries of the Hamiltonian. 
This is described in more detail in Refs. [4] and [5] 
for the gauge glass and XY spin glass, respectively. 
In addition I introduce the link overlap qi defined 

by 

9' = ]J^E(S"-S,")(Sf -Sf), (4) 



with Nf, = dN the number of bonds. 

To avoid critical effects influencing the data, I 
perform the simulations at low temperatures, typ- 
ically T < 0.2Tc. To equilibrate the systems at 
such low temperatures, I use the parallel temper- 
ing Monte Carlo method [7] . For equilibration tests 
for parallel tempering Monte Carlo on the models 
studied I refer the reader to Refs. [4] and [5]. 



Results 

Figures 1 and 2 show data[^for Pdgl) and P{qi) 
at T = 0.050, well below « 0.45 [8] for the 3D 
gauge glass. 




|q| 



Fig. 1. Data for the overlap distribution P(|(?|) at tem- 
perature T = 0.050 for the 3D gauge glass. Note the log- 
arithmic vertical scale. In this and other similar figures in 
the paper, I only display a subset of all the data points 
while the lines connect all the data points in the set. Thus, 
the structure in the lines between neighboring symbols is 
meaningful. 

There is a clear peak in P(|g|) for large \ q\as well 
as a tail at small \q\. The weight in the tail does 
not decrease with increasing L, as is expected in 
the standard interpretations of the droplet theory. 
If anything, the weight increases for larger sizes. 
There is some evidence that f'd?!) stays flat down 
to smaller \q\ for larger L, although the range of 
sizes is too small to make a reliable extrapolation. 



^ I present data for \q\ since for the gauge glass q £ C 
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As with the distribution of \q\, there is a pro- 
nounced peak at large ^/-values in the distribution 
of P{qi) for the 3D gauge glass. Note also the ap- 
pearance of a smaller peak at smaller qi for larger 
system sizes, as predicted by RSB. 



102 



10 - 



0.1 ? 



0.01 







: T = 0.050 




L 




= ■'S 








•5 




6 




1 -8 









0.5 



Fig. 2. The distribution of the link overlap at T = 0.050 
for the 3D gauge glass. Note the logarithmic vertical scale. 



9' + 2{d - ds) = 0.501 ± 0.04. Assuming 6*' « I 
find d-ds= 0.25 ± 0.02, implying that for the 3D 
gauge glass system-size excitations have a fractal 
surface in the thermodynamic limit as predicted 
by the droplet picture. 

Figures 4 and 5 show data for P{q) and P{qi) 
for T = 0.20 (to be compared with Tc « 0.95 [9]) 
for the 4D XY spin glass. Again one sees a large 
peak for large q and qi values. The data for P{qi) 
exhibits a hint of a shoulder for smaller values of qi . 
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Fig. 3. Log-log plot of the variance of qi vs. L at several 
temperatures for the 3D gauge glass. 

The width of the distribution decreases with in- 
creasing system size. This can be seen in Figure 
3 where I plot the variance of the link overlap as 
a function of system size L. The data is consis- 
tent with a power law decrease (a = in Eq. 1) 
where the (presumably effective) exponent varies 
slightly with T. Extrapolating to T = gives fn = 



Fig. 4. Data for the overlap distribution P{q) at T = 0.20 
for 4D XY spin glass. Note the logarithmic vertical scale. 
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Fig. 5. The distribution of the link overlap at T = 0.20 
for the 4D XY spin glass. 

The width of the distribution decreases with in- 
creasing system size. This is demonstrated in Fig. 6 
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where I show the variance of qi against system size 
L for several low temperatures. There is some cur- 
vature in the data for Var((7;) so first I attempt a 
three-parameter fit to Eq. 1 . As there are the same 
number of data points as variables, I cannot assign 
fitting probabilities to the fits but from the data I 
find a = O.OOfOO, 0.00087, 0.00073 and 0.00036 for 
T = 0.200, 0.247, 0.305 and 0.420, respectively. 

I also attempt a power law fit to Eq. f setting 
a = 0. The quality of the fits is poor with proba- 
bilities Q = 5.0 X fO"^ 3.6 X 10-^, 2.9 x 10"^ and 
6.0 X 10-* for T = 0.200, 0.247, 0.305 and 0.42, re- 
spectively]^ The effective exponent is found to 
vary with temperature. Extrapolating to T = 0, I 
obtain^/ ee d' + 2{d-ds) = 0.294 ±0.073. Assum- 
ing 6' « 0, one obtains d ~ ^ O.f 47 ± 0.036. 
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Fig. 6. Log-log plot of the variance of qi vs. L at several 
temperatures for the 4D XY spin glass. 



Conclusions 

To conclude, I have studied the properties of the 
3D gauge glass and the 4D XY spin glass at low 
temperatures. For both models, the order param- 
eter distribution P{q) has, in addition to a peak, 
a tail that seems to extend for smaller values of 
q and whose height seems to persist as the sys- 
tem size increases. This interpretation of the data 



is compatible with the RSB picture or the KMPY 
scenario. The range of lattice sizes is very small, 
however, so this conclusion can be considered at 
most tentative. 

For the gauge glass, the variance of the link 
overlap indicates that the surface of low-energy 
large-scale excitations is fractal with d — dg = 
0.25 ± 0.02 in agreement with the KMPY sce- 
nario. The log- log plot for Yai{qi) for the XY 
spin glass shows curvature, possibly indicating a 
non-zero value in the thermodynamic limit, a re- 
sult compatible with RSB. Due to the small range 
of system sizes, however, I cannot rule out the 
possibility that Yai{qi) ^ at large L, which is 
compatible with the KMPY scenario or droplet 
picture. The effects of vortices in the spin glass 
phase still remain to be understood, ft would be 
useful to look more carefully at the nature of the 
large-scale low-energy excitations to see whether 
they correspond to gradual orientations in the 
spin directions or whether vortices play a role. 
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